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We  introduce,  model  and  solve  to  optimality  a rich  multi-product,  multi-period  and  multi-compartment 
vehicle  routing  problem  with  a required  compartment  cleaning  activity.  This  real-life  application  arises 
in  the  olive  oil  collection  process  in  Tunisia,  where  regional  collection  offices  dispose  of  a fleet  of  vehicles 
to  collect  one  or  several  grades  of  olive  oil  from  a set  of  producers.  For  each  grade,  the  quantity  offered  by 
a producer  changes  dynamically  over  the  planning  horizon.  We  first  provide  a mathematical  formulation 
of  the  problem,  along  with  a set  of  known  and  new  valid  inequalities.  We  then  propose  an  exact  branch- 
and-cut  algorithm  to  solve  the  problem.  We  evaluate  the  performance  of  the  algorithm  on  real  data  sets 
under  different  transportation  scenarios  to  demonstrate  to  our  industrial  partner  the  advantages  of 
using  multi-compartment  vehicles. 
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1.  Introduction 

In  this  paper,  we  introduce,  model  and  solve  a real-world  applica- 
tion of  a multi-product,  multi-period  and  multi-compartment  vehicle 
routing  problem  (MPPC-VRP)  arising  in  the  collection  of  olive  oil  in 
Tunisia.  In  2012,  that  country  was  the  fourth  largest  exporter  of  olive 
oil  worldwide,  with  an  export  production  of  163,000  tons.  This 
amount  was  expected  to  increase  in  2013  according  to  the  General 
Directorate  for  Research  at  the  Ministry  of  Agriculture.  For  climatic  and 
geographical  reasons,  olive  groves  are  rather  widespread  in  the  central 
part  of  the  countiy,  as  shown  in  Fig.  1.  Collecting  olive  oil  is  particularly 
important  during  the  four-month  production  season.  It  mobilizes 
considerable  human  and  material  resources,  and  timeliness  is  crucial 
in  this  operation.  The  producers  work  non-stop  24  h a day  in  order  not 
to  damage  the  harvest.  On  any  given  day,  olive  oil  collection  is  carried 
out  over  six  periods  lasting  almost  24  h in  total.  This  activity  is 
performed  by  a fleet  of  capacitated  heterogeneous  vehicles,  with 
compartments  of  equal  or  different  sizes,  all  equipped  with  a debit 
meter,  enabling  the  decision  maker  to  have  full  knowledge  of  the  load 
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contained  in  each  compartment  at  all  times.  The  oil  must  be  collected 
before  the  producer  runs  out  of  storage  space.  A good  forecast  is 
available  for  the  production  rate  of  each  product  by  each  producer. 

Olive  oil  comes  in  three  different  grades  known  as  extra,  virgin,  and 
lampante.  The  top  two  grades  with  superior  tastes  are  extra  and  virgin, 
which  are  suitable  for  consumption,  whereas  lampante  oil  is  mostly 
destined  for  industrial  uses.  The  transportation  is  regulated  by  law  to 
protect  the  natural  flavors  of  the  oils.  In  particular,  at  each  producer 
site,  a quality  controller  is  in  charge  of  checking  the  oil  grade  proposed. 
Once  the  quality  control  process  has  been  completed  with  success,  the 
quality  controller  seals  the  tank  containing  this  offer.  Thereafter,  once 
the  vehicle  loading  starts,  it  cannot  be  stopped  until  the  tank  is  empty. 
In  addition,  the  different  grades  must  be  kept  separate  during 
transportation,  hence  the  need  to  have  multi-compartment  vehicles. 
It  is  forbidden  to  load  superior  grades  immediately  after  lampante  oil 
in  the  same  compartment,  unless  it  has  been  cleaned  before  the 
changeover.  The  cleaning  activity  incurs  a cost  and  takes  time. 

Routing  problems  with  a cleaning  activity  have  not  been  widely 
studied  from  a scientific  perspective,  but  similar  constraints 
appear  in  other  contexts.  Oppen  et  al.  [33]  consider  the  problem 
of  transporting  different  types  of  live  animals  from  farms  to 
slaughterhouses  by  means  of  multi-compartment  vehicles.  They 
add  time  between  consecutive  tours  to  allow  for  unloading  and 
disinfection  of  the  vehicles.  Hvattum  et  al.  [24]  deal  with  a tank 
allocation  problem  arising  in  the  shipping  of  bulk  oil  and  chemical 
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Fig.  1.  Map  of  Tunisia  pinpointing  producers  and  regional  offices  locations. 
Source:  Google  Maps,  accessed  March  2014. 


products  by  tanker  ships.  They  consider  that  a cleaning  activity  is 
required  if  two  incompatible  products  are  assigned  to  the  same 
compartment  within  less  than  three  trips. 

The  use  of  fleet  with  several  compartments  is  common  in  fuel 
and  oil  distribution  [2,8,10,11,15,38,42]  and  in  some  maritime 
applications  [5,22,24,40],  Transporting  oil  and  fuel  with  multi- 


compartment vehicles  is  more  challenging  and  interesting  from  a 
scientific  point  of  view  than  transporting  food,  where  dry,  refri- 
gerated and  frozen  commodities  can  be  pre-assigned  to  suitable 
compartments.  In  this  case,  the  loading  problem  reduces  to  a 
simple  capacity  checking  procedure  [13,14,30,31].  In  contrast,  in 
fuel  transportation,  a routing  problem  and  a compartment 
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assignment  problem  must  be  solved  jointly.  For  more  details  on 
the  loading  aspect  arising  in  vehicle  routing  problems  with 
compartments,  see  Lahyani  et  al.  [29]  and  Pirkwieser  et  al.  [36], 
Multi-compartment  problems  have  also  been  investigated  in  the 
context  of  inventory-routing  [7,33,35,43],  Other  multi-period 
researches  dealing  with  the  joint  optimization  of  transportation 
and  inventory  include  [1,3,23], 

Multi-compartment  vehicles  are  often  not  equipped  with  debit 
meters,  which  implies  that  a compartment  must  satisfy  at  most  one 
request.  This  scenario  has  been  extensively  studied  [8-11,37,43], 
Coelho  and  Laporte  [6]  introduce  and  define  the  generalized  case 
where  compartments  may  be  equipped  with  debit  meters,  and  the 
load  of  a compartment  may  be  split  between  different  customers. 
They  also  distinguish  between  the  cases  where  customers  may  or 
may  not  be  visited  by  more  than  one  vehicle  per  period.  The 
authors  assess  the  difficulty  of  the  problem  with  split  compart- 
ments and  multiple  visits  per  period.  Derigs  et  al.  [13]  provide  a 
literature  review  on  the  VRP  with  compartments  as  it  arises  in  food 
and  petrol  delivery. 

In  this  paper,  we  introduce  other  complicating  constraints 
reflecting  the  complexities  of  the  olive  oil  industry.  The  mathema- 
tical model  developed  for  this  application  considers  two  vehicle 
types  which  differ  in  the  number  and  sizes  of  the  compartments.  In 
the  olive  oil  industry,  a producer  typically  offers  several  products  in 
different  quantities  in  each  period.  The  grades  and  quantities 
offered  depend  on  the  producer.  All  offers  must  be  picked  up  in 
each  period  since  drivers  are  hired  per  period.  These  assumptions, 
together  with  the  presence  of  the  cleaning  activity,  increase  the 
difficulty  of  the  problem  and  make  it  almost  impossible  to  solve  by 
exact  algorithms  for  most  instances  of  realistic  sizes. 

From  a scientific  perspective,  we  introduce,  model,  and  solve 
exactly  a difficult  and  rich  variation  of  the  well-known  vehicle 
routing  problem  [28],  Part  of  the  complexity  of  the  problem  arises 
from  the  requirement  to  clean  the  compartments,  depending  on 
the  assignment  and  scheduling  decisions.  From  a practical  stand- 
point, we  provide  a tool  that  can  assist  managerial  decision 
making  at  the  tactical  and  operational  levels.  In  particular,  we 
are  able  to  compare  several  transportation  alternatives  for  the 
service  provider  and  to  evaluate  the  potential  routing  savings 
yielded  by  replacing  single  compartment  vehicles  with  multi- 
compartment ones.  Moreover,  we  show  how  our  model  can  be 
extended  to  a number  of  other  industries  operating  under  differ- 
ent assumptions. 

The  remainder  of  this  paper  is  organized  as  follows.  In  Section  2 
we  provide  a formal  description  of  the  problem  and  we  introduce 
a mathematical  model  complemented  with  known  and  new  valid 
inequalities.  The  branch-and-cut  algorithm  is  described  in  Section  3. 
We  present  computational  results  on  real  and  artificially  generated 
instances  in  Section  4,  followed  by  conclusions  in  Section  5. 

2.  Mathematical  description  of  the  problem 

We  first  introduce  some  notation,  followed  by  a mathematical 
model  and  valid  inequalities. 

2.1.  Notation 

The  MPPC-VRP  is  defined  on  an  undirected  graph  Q = (V,  £), 
where  V is  the  vertex  set  and  £ is  the  edge  set.  The  set  V=  (0,  ...,n) 
contains  the  locations  of  the  depot  0 and  of  the  producers 
V'  = {1,  The  set  £ = {(i,j) : ieVj eV',i<j}  defines  the  edges 
of  the  problem.  A routing  cost  ay  is  known  for  each  edge  (ijj  e £. 

The  problem  is  defined  over  a finite  planning  horizon  T = { 1 p}, 

and  a set  K = {1,  ...,K}  of  heterogeneous  vehicles  is  available  at  the 
depot.  Each  vehicle  k contains  a set  Ck  of  compartments  equipped 


with  a debit  meter,  and  each  compartment  l e Ck  has  a capacity  Q,k.  g jS 

The  use  of  vehicle  k incurs  a fixed  cost  (ik.  Identical  products  S ^ 

collected  from  different  producers  can  be  loaded  into  the  same  ' ^ 
compartment  provided  there  is  sufficient  capacity.  Also,  each  ~ £ 
producer  may  be  visited  by  more  than  one  vehicle  in  any  given 
period.  qm 

A set  At  = {1,  ...,M}  of  products  is  offered.  For  ease  of  notation, 
we  assume  that  products  are  ordered  according  to  their  grade,  and  if 
the  product  of  the  lowest  grade,  i.e„  M,  is  loaded  into  a compartment 
of  a vehicle,  it  contaminates  the  compartment  and  a cleaning 
procedure  is  required  to  load  any  product  m<  M (me  M)  in  sub- 
sequent periods.  If  a cleaning  procedure  is  undertaken,  it  generates  a 
cost  y.  Each  producer  i e V offers  a quantity  d[m  of  product  of  type  m 
in  period  t with  which  is  associated  a transportation  request.  Not  all 
producers  offer  all  types  of  products,  and  the  available  quantities 
change  from  one  period  to  the  next. 

2.2.  Mathematical  model 

We  now  provide  an  integer  linear  programming  formulation 
for  the  MPPC-VRP,  followed  by  known  and  new  valid  inequalities 
used  to  strengthen  the  mathematical  model.  The  decision  vari- 
ables are  defined  as  follows.  The  routing  variables  xjf  are  equal  to 
the  number  of  times  edge  (ij)  is  used  by  vehicle  k in  period  t;  the 
visiting  variables  ykt  are  binary  and  equal  to  one  if  and  only  if 
vertex  i is  visited  by  vehicle  k in  period  t;  the  assignment  variables 
zjm  are  equal  to  one  if  and  only  if  product  m of  producer  i is  loaded 
into  compartment  / of  vehicle  k in  period  t;  likewise,  the  variables 
wj'  are  equal  to  one  only  if  product  m is  loaded  in  compartment 
/ of  vehicle  k in  period  t.  Finally,  the  variables  ulkt  are  binary  and 
indicate  whether  or  not  compartment  / of  vehicle  k is  cleaned  at 
the  beginning  of  period  t. 

Recall  that  not  all  producers  offer  all  products  in  each  period. 
Nevertheless,  we  define  variables  zjfjj  for  all  combinations  of 
producers  and  products,  because  this  simplifies  the  formulation. 

If  a producer  i does  not  offer  product  m in  period  t,  we  then  set 
zjjjj  = 0 for  all  compartments  / and  vehicles  k.  Finally,  we  need  to 
define  the  set  T\  containing  the  periods  elapsed  between  s and  t, 
i.e„  T\  = {s+l. ....  t- 1). 

The  problem  can  then  be  modeled  as  follows: 

minimize  2 I 2 I 2 I ZPVo+Z  I 2 yulkt 


teTke  K(ij)  £ £ tsTkeK jeV  tsTkeKlEck 

(1) 

subject  to 

2 <4,zjm  < Q."1  keK,  leCk,  meM,  teT  (2) 

ieV 

z^<ykt  ieV' , keK.,  leCk,  meM,  teT  (3) 

z!m^wrn  ieV,  keK,  leCk,  meM,  teT  (4) 

wlkt  < 2 z!m  keK,  leCk,  meM,  teT  (5) 

ieV 

2 w“<l  keK,  l eCk,  teT  (6) 

meM 

2 xf+  2 *jf  = 2yf  ieV,  keK,  teT  (7) 

j e V,i  < j j e Vj  < i 

2 2 xjf<  2 ZsV,  ZeZ,  keK,  teT  (8) 

i e Zj  e Z,i  <j  i e Z 

2 2z!m  = 1 ieV',  m e M,  teT  (9) 

ksKlgrt 
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wlM+wlm<\+ulkt+  2 ulkt'  ke/C,  leCk , meM\[M}, 

t'eT's 

S,teT  (10) 

Xgj  e {0, 1,2}  jeV,  fee/c,  teT  (11) 

Xjf  e {0, 1 } ij  e V',  he  1C,  teT  (12) 

ykt,Zj£,  wjjf,  ulkt  e {0, 1)  ieV,  jeV,  kel C,  1 e £fc,  meM,  teT. 

(13) 

The  objective  function  (1)  minimizes  the  sum  of  the  routing 
cost,  the  vehicle  fixed  cost  and  the  compartment  cleaning  cost. 
Constraints  (2)  ensure  that  the  capacity  of  each  compartment  is 
not  violated.  Constraints  (3)— (5)  link  the  variables  ykt,  and  zj^. 

Specifically,  constraints  (3)  ensure  that  a product  from  a producer 
is  loaded  into  a given  compartment  of  a vehicle  at  a given  period 
only  if  the  producer  is  served  by  the  vehicle.  Constraints  (4)  and 
(5)  guarantee  that  a compartment  is  allowed  to  carry  a product  in 
a given  period  only  if  the  vehicle  visits  a producer  offering  that 
product  in  the  same  period.  Constraints  (6)  ensure  that  each 
compartment  carries  at  most  one  type  of  product  at  a time. 
Constraints  (7)  and  (8)  are  degree  and  subtour  elimination  con- 
straints, respectively.  Constraints  (9)  ensure  that  all  the  quantities 
being  offered  are  collected.  Constraints  (10)  mean  that  a cleaning 
operation  is  performed  if  necessary.  The  term  £t,  sTtU!kt  keeps 
track  of  each  compartment  cleaning  operations  for  the  interval 
between  periods  s and  t.  Constraints  (11 )— (13)  define  the  integer 
and  binary  nature  of  the  variables. 

2.3.  Model  extension 

The  model  defined  by  ( 1 )— ( 13)  can  be  generalized  and 
extended  to  formulate  other  problems  or  to  handle  different 
assumptions.  For  example,  this  formulation  is  valid  for  the  food 
distribution  problem  [4,13]  or  for  waste  collection  [31]  where 
different  commodities  can  be  pre-assigned  to  suitable  compart- 
ments. In  many  other  cases,  small  ad  hoc  modifications  suffice  to 
model  the  problem  using  a very  similar  formulation.  In  the  current 
setting,  we  assume  that  different  requests  for  the  same  product 
can  be  loaded  into  the  same  compartment  and  that  the  producers 
may  receive  more  than  one  visit  by  different  vehicles  in  any 
period.  These  assumptions  characterize  the  general  case  which 
includes  an  extra  layer  of  difficulty  (brought  by  the  cleaning 
activities)  for  the  multi-compartment  vehicle  routing  problem 
(MC-VRP).  However,  vehicles  are  not  always  equipped  with  debit 
meters,  as  is  the  case  in  the  fuel  industry  [8-11  ].  In  other  cases,  the 
producer  may  impose  a limit  of  at  most  one  visit  per  time  period 
[8,39,41],  Coelho  and  Laporte  [6]  showed  how  to  handle  such 
specific  versions  of  the  MC-VRP. 

In  order  to  account  for  these  changes  in  our  current  model,  no 
new  decision  variable  is  needed.  We  can  simply  remove  the  sum 
over  i in  constraints  (2)  to  respect  compartment  capacity  and  add 
some  side  constraints  similar  to  (6).  For  example,  constraints  (14) 
limit  the  maximum  number  of  visits  for  a given  producer,  and 
constraints  (15)  prevent  the  load  of  a compartment  from  being 
split: 

1 ieV,  teT  (14) 

ke/C 

2 2z'm<l  ke/C.,  leC,  teT.  (15) 

me  Mie  V' 

Another  assumption  of  the  problem  at  hand  is  that  we  have  no 
control  over  the  inventory  levels  at  producers  and  all  producers 
have  the  same  priority.  However,  in  other  applications  such  as  the 


petrol  station  replenishment  problem  [10]  and  the  livestock 
collection  problem  [32,33],  the  decision  maker  may  have  different 
priorities  for  each  producer.  To  this  end,  at  the  beginning  of  each 
time  period,  an  extra  set  of  variables  would  identify  producers 
with  a critical  storage  level  and  set  their  priority  higher,  which 
would  then  be  reflected  in  the  objective  function  in  order  to  create 
a prioritization  of  the  visits. 

The  problem  we  have  modeled  assumes  that  each  request  is 
fully  assigned  to  one  compartment.  In  some  applications, 
splitting  a request  and  loading  it  into  different  compartments  or 
vehicles  is  allowed  [6,13,32,33].  In  such  cases,  variables  z[m  are 
continuous  and  represent  the  fraction  of  product  m associated 
with  producer  i being  loaded  in  compartment  / of  vehicle  k.  This 
extension  requires  modifications  to  constraints  (5).  Constraints  (2)- 
(4)  and  (9)  remain  valid. 

Another  variant  of  this  problem  is  that  arising  in  the  retail 
delivery  sector,  where  the  items  loaded  into  compartments  are 
measured  not  only  by  their  volume,  but  by  their  precise  dimen- 
sions and  weight  [16,17,20,25],  In  our  existing  model,  new  con- 
straints similar  to  (2)  would  be  required  to  handle  this  case. 

Finally,  we  have  modeled  the  problem  assuming  the  vehicle 
fleet  is  heterogeneous,  and  that  each  vehicle  contains  an  hetero- 
geneous set  of  compartments.  Obviously,  the  homogeneous  ver- 
sion of  the  problem  is  a special  case  of  ours  [6,13,14],  and  our 
formulation  remains  valid  for  it. 

2.4.  Valid  inequalities 

The  formulation  defined  by  (1 )— (13)  is  sufficient  to  model  the 
MPPC-VRP.  We  can,  however,  strengthen  it  through  the  inclusion 
of  valid  inequalities  in  the  form  of  symmetry  breaking  constraints 
and  additional  cuts  imposing  bounds  on  the  integer  variables.  The 
first  ones  are  related  to  the  period  in  which  a cleaning  operation 
takes  place.  For  example,  suppose  a contaminating  product  is 
transported  by  a given  compartment  in  period  t=  1,  and  the  next 
use  of  this  compartment  is  to  carry  a higher  grade  product  in 
period  t=  5.  Then,  at  least  four  optimal  solutions  exist,  by  cleaning 
the  compartment  in  period  2,  3,  4,  or  5.  In  order  to  avoid  such 
symmetries,  we  impose  constraints  (16)  which  postpone  the 
cleaning  operation  as  much  as  possible: 

Ulkt  < 2 wjjf  ke  K,le  Ck,t  eT.  (16) 

meM 

We  also  integrate  the  vehicle  and  the  compartment  symmetry 
breaking  constraints  for  the  first  period,  when  no  cleaning  opera- 
tion is  necessary.  These  constraints  are  inspired  from  those 
proposed  in  Coelho  and  Laporte  [6],  Constraints  (17)  and  (18)  are 
valid  if  the  considered  vehicles  and  compartments  are  identical. 
We  define  the  set  K!  c K.  containing  only  the  homogeneous 
vehicles  of  K..  Note  that  in  the  olive  oil  industry  as  in  petrol 
distribution,  vehicles  with  different  compartments  capacities  are 
seldom  used  in  order  to  reduce  the  imbalance  of  loaded  vehicles 
on  the  road  [13].  These  symmetry  breaking  constraints  are 

ZTi's  lyf-1'1  keK?\{  1)  (17) 

i e V i e V 

and 

w“  < 2 ke/C,le£k\{lj,heM.  (18) 

meM 

Constraints  (17)  rank  identical  vehicles  according  to  the  index  of 
the  producers  served.  In  particular,  they  ensure  that  among 
vehicles  of  the  same  type,  vehicle  k cannot  serve  any  customer  if 
vehicle  fc-1  has  not  already  been  used  in  the  first  period. 
Constraints  (18)  rank  identical  compartments  of  vehicle  k.  They 
state  that  if  deliveries  are  performed  using  compartment  /,  then 
compartment  /-I  has  already  been  used.  These  rules  cannot  be 
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generalized  to  the  remaining  periods  because  they  may  impact  the 
cleaning  operation  and  ultimately  affect  the  solution  cost. 

We  also  make  use  of  a known  set  of  useful  cuts  to  enforce 
logical  relationships  between  routing,  visiting  and  assignments 
variables.  For  more  details  on  logical  inequalities  for  routing 
problems,  see  Coelho  and  Laporte  [6]  and  Gendreau  et  al.  [19], 
These  cuts  are  as  follows: 


xoi  ^ 2y  f i eV',  ke/C,  teT 

(19) 

Xjf  < yfc  ij  eV,  ke/C,  teT 

(20) 

yfC  ^yo  ieV,  k e K,  teT 

(21) 

yf<  Z Z45  ieV,  ke/c , 

meMleCt 

teT 

(22) 

y?  < Z Z wm  ' e v'>  k e 

m e Mj  E £*< 

teT 

(23) 

ZfNZ  Z Z4£  ke/c. 

teT. 

(24) 

i e V ie  V'm  e Ml  e £* 


Constraints  (19),  (20)  and  (21 ) are  referred  to  as  routing  cuts,  as 
a way  to  ensure  that  if  a producer  i is  visited,  i.e.,  (19)  holds,  or  a 
producer  j is  the  successor  of  a producer  i,  i.e.,  (20)  holds  on  the 
route  of  vehicle  k in  period  t,  then  producer  i must  be  visited,  i.e., 
y\c  = 1.  Constraints  (20)  can  also  be  considered  as  subtour  elim- 
ination constraints  (8)  when  \Z\  = 2.  Inequalities  (21)  guarantee 
that  if  vehicle  k visits  producer  i in  period  f,  then  the  depot  must 
be  included  in  the  route  of  vehicle  k in  period  t.  Through 
constraints  (22)  and  (23),  we  ensure  that  if  a producer  i is  visited 
in  period  t by  vehicle  k,  then  at  least  one  product  m should  be 
loaded  in  some  compartment  of  vehicle  k in  that  period.  Con- 
straints (24)  strengthen  the  relationships  between  the  collection 
routes,  products  and  compartments.  Specifically,  a collection  route 
using  vehicle  k in  period  t exists  to  ensure  the  pickup  of  some 
products  from  a producer  i and  load  them  in  some  compartment 
of  that  vehicle.  We  also  note  that  constraints  (24)  are  the  sum  over 
the  locations  of  (22).  Even  if  these  constraints  are  redundant  in 
this  context,  they  are  known  to  help  the  mathematical  program- 
ming solvers  derive  new  cuts  and  improve  the  overall  algorithmic 
performance  [7,21,26,27], 


3.  Branch-and-cut  algorithm 

We  have  implemented  a branch-and-cut  algorithm  capable  of 
solving  the  formulation  just  presented.  All  variables  of  the  for- 
mulation are  explicitly  handled  by  the  algorithm.  Since  the 
number  of  constraints  (16)-(24)  is  polynomial,  they  are  added  a 
priori  to  the  model.  In  the  sequel,  we  will  show  how  each  subset  of 
inequalities  impacts  its  solution.  In  contrast,  we  cannot  generate 
all  subtour  elimination  constraints  (8)  a priori  since  their  number 
is  exponential.  These  are  dynamically  generated  as  cuts  as  they  are 
found  to  be  violated.  This  procedure  is  described  in  Section  3.1. 
The  formulation  is  then  solved  by  branch-and-cut  as  follows.  At  a 
generic  node  of  the  search  tree,  a linear  program  with  relaxed 
integrality  constraints  is  solved,  a search  for  violated  constraints  is 
performed,  appropriate  valid  inequalities  are  added  to  eliminate 
subtours,  and  the  current  subproblem  is  then  reoptimized.  This 
process  is  reiterated  until  a feasible  or  dominated  solution  has 
been  reached,  or  until  no  more  cuts  can  be  added.  At  this  point, 
branching  on  a fractional  variable  occurs.  We  provide  in  Algorithm  1 a 
sketch  of  the  branch-and-cut  scheme. 


ft 

r*  0 

3.1.  Separation  algorithm  for  subtour  elimination  constraints 

« I 

When  the  linear  program  without  subtour  elimination  con-  ? ^ 
straints  (8)  is  solved  and  an  integer  solution  is  obtained,  two  ” i 
situations  can  occur.  The  first  one  consists  of  a solution  in  which  cgl 
no  subtours  are  present,  and  the  solution  is  then  feasible.  In  the 
second  case,  one  can  attempt  to  identify  the  number  of  connected 
components  by  means  of  different  algorithms,  namely  by  comput- 
ing the  minimum  cut  in  the  graph  or  its  equivalent  problem,  the 
maximum  flow  [12],  The  procedure  described  in  Padberg  and 
Rinaldi  [34]  consists  of  constructing  an  auxiliary  graph  as  follows. 

First,  all  nodes  visited  in  the  original  problem  are  added  to  the 
auxiliary  graph.  Then,  the  value  of  each  routing  variable  associated 
with  the  partial  solution  is  added  as  a capacity  on  the  edge  linking 
the  two  corresponding  nodes  in  the  auxiliary  graph.  Finally,  a 
max-flow  problem  between  the  depot  and  each  node  is  solved.  If 
the  max-flow  is  equal  to  zero,  this  means  that  the  node  belongs  to 
a component  which  is  not  connected  to  the  depot.  By  doing  this, 
one  can  already  identify  two  sets,  S being  the  set  containing  the 
depot,  and  its  complement  S with  all  the  remaining  nodes  not 
connected  to  the  depot.  Alternatively,  by  selecting  different  origins 
instead  of  the  depot  and  running  a max-flow  algorithm  for  each 
other  node,  one  can  easily  identify  all  the  connected  components 
of  this  auxiliary  graph.  Once  we  have  identified  that  more  than 
one  connected  component  exists,  we  need  to  add  cuts  to  forbid 
such  a solution.  This  is  achieved  by  adding  the  appropriate 
constraints  (8)  related  to  the  nodes  of  S. 

This  cutting  procedure  is  embedded  within  a branch-and- 
bound  scheme.  At  each  node  of  the  branch-and-bound  tree  we 
evaluate  the  current  solution  for  subtours,  separate  the  appro- 
priate constraints,  add  them  to  the  model,  and  reoptimize  the 
node  until  no  more  violated  subtour  elimination  constraints  can 
be  identified. 

Algorithm  1.  Branch-and-cut  algorithm. 

1 : At  the  root  node  of  the  search  tree,  generate  and  insert  all 

valid  inequalities  into  the  program. 

2:  z*<-oo. 

3:  Termination  check: 

4:  if  there  are  no  more  nodes  to  evaluate  then 

5:  Stop  with  the  incumbent  and  optimal  solution  of  cost  z*. 

6:  else 

7:  Select  one  node  from  the  branch-and-bound  tree. 

8:  end  if 

9:  Subproblem  solution:  solve  the  LP  relaxation  of  the  node 

and  let  z be  its  cost. 

10:  if  the  current  solution  is  feasible  then 

11:  if  z > z*  then 

12:  Go  to  the  termination  check. 

13:  else 

14:  z*<-z. 

15:  Update  the  incumbent  solution. 

16:  Prune  the  nodes  with  a lower  bound  larger  than  or 

equal  to  z*. 

17:  Go  to  the  termination  check. 

18:  end  if 

19:  end  if 

20:  Cut  generation: 

21 : if  the  solution  of  the  current  LP  relaxation  violates  any  cuts 

then 

22:  Identify  connected  components  as  in  Padberg  and 

Rinaldi  [34], 

23:  Determine  whether  the  component  containing  the 

producer  is  weakly  connected  as  in  Gendreau  et  al.  [18], 

24:  Add  violated  subtour  elimination  constraints  (8). 
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25:  Go  to  the  subproblem  solution. 

26:  end  if 

27:  Branching:  branch  on  one  of  the  fractional  variables. 

28:  Go  to  the  termination  check. 

4.  Computational  experiments 

In  this  section  we  describe  the  computational  experiments 
carried  out  in  order  to  assess  the  performance  of  our  model  and 
algorithm.  We  provide  in  Section  4.1  details  of  the  real  instances  we 
have  obtained  from  our  Tunisian  partner,  and  the  instances  gener- 
ated with  a different  fleet  composition.  In  Section  4.2  we  describe 
the  results  of  computational  experiments  performed  to  evaluate  the 
effectiveness  of  the  cuts  and  valid  inequalities,  and  we  compare  our 
solutions  with  those  corresponding  to  the  current  situation. 

We  have  coded  the  algorithm  in  C+  + using  IBM  CPLEX  Studio 
12.5.1  as  the  MIP  solver.  All  computations  were  executed  on  a grid 
of  Intel  Xeon™processors  running  at  2.66  GHz  with  up  to  48  GB  of 
RAM  installed  per  node,  with  the  Scientific  Linux  6.1  operating 
system.  All  instances  and  detailed  results  are  available  from  http:// 
www.leandro-coelho.com. 

4.1.  Instance  generation 

We  have  created  a set  of  five  instances  based  on  real  data 
gathered  from  industrial  partners  in  the  regions  of  Sfax  and 
Kairouan  in  Tunisia.  In  terms  of  size,  we  handle  instances  with 
one  depot  and  up  to  45  transportation  requests  loaded  in  three  or 
four  vehicles.  The  product  quantities  being  offered  are  either 


Table  1 

Producers  locations  in  Figs.  2 and  3. 


Producer 

Location 

Producer 

Location 

1 

Feriana,  Kasserine 

A 

Bir  El  Hfey,  Sidi  Bouzid 

2 

Jelma,  Sidi  Bouzid 

B 

Sened,  Gafsa 

3 

Regueb,  Sidi  Bouzid 

C 

Sidi  Ali  Ben  Aoun,  Sidi  Bouzid 

4 

Hajeb  El  Ayoun,  Kairouan 

D 

North  of  Sfax 

5 

El  Khazaziya,  Kairouan 

E 

Road  Gremda,  Sfax 

6 

Boussari,  Kairouan 

F 

Road  Mahdia,  Sfax 

7 

Cherarda,  Kairouan 

G 

Sidi  Bouzid 

8 

El  Houareb,  Kairouan 

9 

South  of  Kairouan 

obtained  from  our  partner,  or  estimated  when  these  could  not 
be  made  publicly  available  for  confidentiality  reasons.  There  is  no 
restriction  on  the  number  of  producers  that  can  be  visited  on  a 
route.  The  quantities  of  available  extra,  virgin  and  lampante  oils 
represent  56%,  30%  and  14%  of  the  total  available  production, 
respectively.  The  routable  network  dataset  is  constructed  using 
real  travel  distances.  Only  driving  distances  provided  by  the 
regional  office  of  Sfax  take  traffic  congestion  and  the  road  state 
into  account.  The  geographical  locations  of  the  producers  of  this 
case  study  are  listed  in  Table  1 and  depicted  in  Figs.  2 and  3. 

The  industrial  partners  of  the  regions  Sfax  and  Kairouan  do  not 
currently  dispose  of  the  same  fleet  composition.  The  first  office 
controls  a heterogeneous  limited  size  fleet  of  single-compartment 
vehicles  having  a capacity  of  10  tons,  and  double-compartment 
vehicles  in  which  each  compartment  has  a capacity  of  five  tons. 
These  two  configurations  will  be  denoted  by  type  1 and  type  II, 
respectively.  The  regional  office  of  Kairouan,  which  presently  uses 
single-compartment  vehicles  only,  will  revise  its  short-term  pro- 
curement policy  after  we  will  have  identified  the  cost  savings 
achieved  by  using  multi-compartment  vehicles.  Some  drivers  are 
outsourced,  which  enables  the  service  providers  to  perform  tours 
and  to  exploit  the  available  fleet  over  six  periods  spread  out  during 
day  time  and  in  the  evening.  The  data  exploitation  results  in  five 
original  instances,  with  an  asterisk  prepended  to  their  names.  To 
cover  the  different  scenarios  under  a different  fleet  composition, 
we  generate  10  instances  identical  to  the  original  ones,  but  with  a 
different  set  of  vehicles.  The  names  of  the  test  instances  highlight 
the  factors  that  may  affect  the  results.  These  factors  include  the 
first  four  letters  of  the  regional  office  they  refer  to,  the  total 
number  of  requests,  and  the  number  of  vehicles  of  type  I and  type  II, 
respectively.  Since  the  original  instances  may  contain  some  pro- 
ducts with  a supply  of  more  than  five  tons,  some  restrictions  had  to 
be  made  while  generating  the  fleet  composition.  Each  collection  route 
requires  at  least  one  single-compartment  vehicle  for  the  Kairouan 
region,  and  two  single-compartment  vehicles  for  the  Sfax  region. 
Table  2 summarizes  the  characteristics  of  these  instances.  Regarding 
the  objective  function,  we  have  used  the  following  parameters  after  a 
tuning  phase  and  discussions  with  our  industrial  partner: 

• a,.  = € I per  driven  kilometer  between  vertices  i and  j; 

• P = € 15  per  vehicle  k used  per  period  t regardless  of  the 
vehicle  type; 

• ^ = €10  per  cleaning  activity. 
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Fig.  2.  Geographical  locations  of  the  producers  around  the  region  of  Sfax. 
Source : Google  Maps,  accessed  March  2014). 
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Fig.  3.  Geographical  locations  of  the  producers  around  the  region  of  Kairouan. 
Source:  Google  Maps,  accessed  March  2014. 


Table  2 

Configurations  of  the  real-world  instances. 


Instance 

# Producers 

# Requests 

Fleet  composition 
Type  I Type  II 

*I<air  27  3 0 

6 

27 

3 

0 

Kair_27_l_2 

6 

27 

1 

2 

Kair  27  2 1 

6 

27 

2 

1 

*I<air_33_3_0 

5 

33 

3 

0 

Kair_33_l_2 

5 

33 

1 

2 

Kair  33  2 1 

5 

33 

2 

1 

*I<air  34  3 0 

6 

34 

3 

0 

Kair_34_l_2 

6 

34 

1 

2 

Kair  34  2 1 

6 

34 

2 

1 

*Kair_45_3_0 

7 

45 

3 

0 

Kair_45_l_2 

7 

45 

1 

2 

Kair  45  2 1 

7 

45 

2 

1 

*Sfax  39  2 2 

7 

39 

2 

2 

Sfax  39  3 1 

7 

39 

3 

1 

Sfax_39_4_0 

7 

39 

4 

0 

4.2.  Computational  results 

We  have  run  the  proposed  algorithm  over  the  data  sets  shown 
in  Table  2.  Table  3 summarizes  the  performance  of  the  algorithm 
on  the  15  real  instances.  We  assess  the  performance  of  the 
symmetry  breaking  cuts  and  the  routing  and  assignment  cuts  by 
comparing  the  solutions  obtained  for  different  configurations.  The 
implementation  with  constraints  (1 )— (13)  is  used  as  a reference 
point.  Specifically,  for  each  instance  we  present  the  number  of 
nodes  explored  in  the  branch-and-cut  tree,  the  ratio  of  the  lower 
bound  at  the  root  node  between  the  configuration  with  cuts  and 
the  basic  configuration  (1 )— (13),  and  finally  the  running  time  in 
seconds. 

The  algorithm  proves  optimality  over  all  the  15  instances 
within  very  short  computing  times.  For  most  of  the  instances  of 
the  region  of  Kairouan  with  three  vehicles  and  up  to  seven 


producers,  the  algorithm  takes  less  than  one  second  to  reach 
optimality.  However,  it  requires  more  computational  effort  for  the 
instances  of  Sfax,  especially  when  the  proposed  cuts  are  disabled. 
This  may  be  explained  by  the  fact  that  increasing  the  number  of 
vehicles  generates  much  more  symmetry. 

A deeper  analysis  of  the  tested  configurations  shows  that  the 
introduction  of  valid  inequalities  significantly  improves  the  per- 
formance of  the  algorithm.  The  average  running  time  is  reduced 
from  36.53  to  1.46  s and  the  average  number  of  explored  nodes 
goes  down  from  8143  to  211  when  the  full  model  is  implemented. 
The  short  computational  time  results  from  the  fact  that  the  model 
provides  a high  quality  lower  bound  at  the  root  node  of  the  search 
tree.  On  average,  the  lower  bound  of  the  model  with  all  cuts  is 
almost  equal  to  1.5  times  the  initial  lower  bound  value  of  the  basic 
formulation  (1 )— (13).  The  full  set  of  proposed  cuts  is  essential  to 
achieve  the  best  algorithmic  performance.  In  particular,  comparing 
the  first  two  configurations  substantiates  the  efficiency  of  the 
symmetry  breaking  constraints  (16)-(18).  Even  though  the  quality 
of  the  lower  bound  is  not  improved,  the  introduction  of  these 
constraints  reduces  the  number  of  explored  nodes  in  the  branch- 
and-cut  tree  on  average  from  8143  to  5685  and  reduces  the 
running  time  for  most  of  the  instances.  The  assignment  cuts 
(22)-(24)  have  a more  positive  impact  on  these  instances  than 
the  routing  cuts  (19)-(21).  Without  the  assignment  cuts,  the 
model  requires  more  iterations  and  explores  more  nodes  to  prove 
optimality.  On  average,  the  algorithm  explores  respectively  823 
and  3740  nodes  when  using  separately  assignment  cuts  and 
routing  cuts.  Under  these  two  configurations,  we  have  obtained 
almost  the  same  improvement  in  the  lower  bound  value  (1.31 
against  1.20)  by  imposing  the  assignment  cuts. 

When  this  study  was  undertaken,  our  industrial  partner  was 
trying  to  find  ways  of  minimizing  its  logistics  costs,  i.e.,  the  fixed 
and  variable  costs  related  to  transportation,  as  well  as  the  related 
cleaning  costs.  In  particular,  it  was  paying  close  attention  to  the 
fleet  composition  component  and  its  impact  on  the  overall  costs. 
Using  our  methodology,  we  could  provide  alternative  solutions  to  the 
managers  by  generating  10  instances  while  the  overall  capacity 
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Table  3 

Summary  of  computational  results  of  the  test  instances  with  respect  to  several  configurations. 


Model 

Instance 

# Nodes 

Lower  bound  increase 

Time(s) 

(1)-(13) 

*Kair_27_3_0 

3674 

1 

9 

(1 )— (13),  (16)— (18) 

1051 

1 

3 

(1)— (13),  (19)— (21) 

171 

1.28 

1 

(1 )— (13),  (22)-(24) 

101 

1.45 

1 

(1 )— (13),  (16)— (IS),  (19)-(21),  (22)-(24) 

2 

1.54 

0 

(1)— (13) 

Kair_27_l_2 

620 

1 

7 

(1 )— (13),  (16)— (18) 

1348 

1 

11 

(1)— (13),  (19)— (21) 

227 

1.28 

2 

(1 )— (13),  (22)— (24) 

164 

1.13 

2 

(1 )— (13),  (16)— (IS),  ( 19)— (21 ),  (22)-(24) 

46 

1.54 

1 

(1)— (13) 

Kair_27_2_l 

96 

1 

3 

(1 )— (13),  (16)— (IS) 

54 

1 

2 

(1 )— (13),  (19)— (21 ) 

41 

1.25 

2 

(1 )— (13),  (22)— (24) 

11 

1.29 

1 

(1 )— (13),  (16)— (IS),  ( 19)— (21 ),  (22)-(24) 

1 

1.39 

0 

(1)— (13) 

*Kair_33_3_0 

190 

1 

1 

(1 )— (13),  (16)— (IS) 

274 

1 

1 

(1 )— (13),  (19)— (21 ) 

5 

1.57 

0 

(1 )— (13),  (22)-(24) 

14 

1.36 

1 

(1 )— (13),  (16)— (18),(19)— (21 ),  (22)-(24) 

1 

1.80 

0 

(1 )— (13) 

Kair_33_l_2 

141 

1 

2 

(1 )— (13),  (16)— (IS) 

14 

1 

2 

(1)— (13),  (19)— (21) 

44 

1.40 

1 

(1)— (13),  (22)-(24) 

95 

1.13 

1 

(1 )— (13),  (16)-(18),(19)-(21),  (22)-(24) 

84 

1.43 

0 

(1)— (13) 

Kair_33_2_l 

44 

1 

1 

(1 )— (13),  (16)— (18) 

42 

1 

2 

(1)— (13),  (19)— (21) 

48 

1.51 

1 

(1 )— (13),  (22)— (24) 

30 

1.25 

1 

(1 )— (13),  (16)— (18),(19)— (21),  (22)— (24) 

15 

1.68 

1 

(1)— (13) 

*Kair_34_3_0 

1968 

1 

6 

(1 )— (13),  (16)— (IS) 

915 

1 

3 

(1)— (13),  (19)— (21 ) 

231 

1.28 

2 

(1 )— (13),  (22)-(24) 

61 

1.43 

2 

(1 )— (13),  (16)— (1S),(19)— (21),  (22)— (24) 

2 

1.53 

0 

(1)— (13) 

Kair_34_l_2 

205 

1 

2 

(1)— (13),  (16)— (IS) 

186 

1 

2 

(1)— (13),  (19)— (21 ) 

74 

1.28 

1 

(1 )— (13),  (22)-(24) 

204 

1.13 

2 

(1 )— (13),  (16)— (18),  (19)-(21),  (22)-(24) 

27 

1.29 

0 

(1)— (13) 

Kair_34_2_l 

345 

1 

3 

(1)— (13),  (16)— (18) 

564 

1 

3 

(1 )— (13),  (19)— (21 ) 

134 

1.25 

2 

(1 )— (13),  (22)— (24) 

76 

1.29 

2 

(1 )— (13),  (16)— (1S),(19)— (21 ),  (22)-(24) 

20 

1.40 

1 

(1)— (13) 

*Kair_45_3_0 

653 

1 

5 

(1 )— (13),  (16)— (IS) 

884 

1 

6 

(1 )— (13),  (19)— (21 ) 

165 

1.20 

3 

(1 )— (13),  (22)-(24) 

12 

1.18 

1 

(1 )— (13),  (16)— (1S),(19)— (21),  (22)— (24) 

1 

1.39 

0 

(1)— (13) 

Kair_45_l_2 

382 

1 

6 

(1 )— (13),  (16)— (IS) 

151 

1 

4 

(1 )— (13),  (19)— (21 ) 

176 

1.21 

3 

(1 )— (13),  (22)-(24) 

95 

1.08 

4 

(1)-(13),  (16)— (1S),(19)— (21 ),  (22)— (24) 

180 

1.26 

2 

(1)— (13) 

Kair_45_2_l 

446 

1 

6 

(1 )— (13),  (16)— (IS) 

229 

1 

5 

(1)— (13),  (19)— (21 ) 

33 

1.19 

2 

(1 )— (13),  (22)-(24) 

38 

1.12 

1 

(1 )— (13),  (16)-(18),(19)-(21),  (22)-(24) 

28 

1.31 

1 

(1)— (13) 

*Sfax  39  2 2 

24,736 

1 

119 

(1 )— (13),  (16)— (IS) 

17,390 

1 

83 

(1)— (13),  (19)— (21) 

851 

1.37 

6 

(1)— (13),  (22)-(24) 

1920 

1.10 

10 

(1 )— (13),  (16)— (18),(19)— (21 ),  (22)-(24) 

415 

1.85 

3 

(1)— (13) 

Sfax_39_3_l 

19,888 

1 

92 

(1 )— (13),  (16)— (IS) 

17,692 

1 

79 

(1)— (13),  (19)— (21) 

18,726 

1.29 

75 

(1 )— (13),  (22)— (24) 

5395 

1.04 

24 

(1 )— (13),  (16)— (1S),(19)— (21 ),  (22)— (24) 

1301 

1.77 

7 

(1)— (13) 

Sfax  39  4 0 

68,780 

1 

285 

(1 )— (13),  (16)— (IS) 

44,506 

1 

181 

(1 )— (13),  (19)— (21 ) 

35,202 

1.31 

123 

(1 )— (13),  (22)-(24) 

4145 

1.04 

17 

(1)— (13),  (16)— (1S),(19)— (21),  (22)-(24) 

1071 

1.79 

6 

ifV 
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5.  Conclusions 


.36 
2 I 


Table  4 

Relative  savings  compared  to  existing  situation. 


Instance 

Optimized  solution 

Manual  solution 

# Km 

# Vehicles 

# Cleanings 

# Km 

# Vehicles 

# Cleanings 

*Kair  27  3 0 

2361 

17 

1 

2435 

17 

0 

Kair_27_l_2 

1848 

15 

0 

- 

- 

- 

Kair_27_2_l 

1848 

15 

1 

- 

- 

- 

*I<air  33  3 0 

2173 

16 

0 

2339 

16 

0 

Kair_33_l_2 

1990 

16 

0 

- 

- 

- 

Kair_33_2_l 

2104 

16 

0 

- 

- 

- 

*I<air  34  3 0 

2737 

17 

1 

2829 

17 

1 

Kair_34_l_2 

2322 

16 

1 

- 

- 

- 

Kair_34_2_l 

2424 

17 

1 

- 

- 

- 

*Kair  45  3 0 

4477 

18 

0 

4721 

18 

0 

Kair_45_l_2 

4008 

17 

3 

- 

- 

- 

Kair  45  2 1 

4388 

17 

0 

- 

- 

- 

*Sfax_39_2_2 

2755 

21 

0 

2866 

21 

0 

Sfax  39  3 1 

2787 

21 

0 

- 

- 

- 

Sfax_39_4_0 

2811 

21 

0 

- 

- 

- 

remains  unchanged  and  the  number  of  vehicles  of  types  I and  11  varies. 
Table  4 summarizes  these  results.  We  note  that  substantial  savings  are 
achieved  if  both  types  of  vehicles  are  used.  We  observe  that  a 
combination  of  vehicles  yields  better  quality  solutions  for  all  Kairouan 
instances,  with  improvements  ranging  from  1.9%  to  21.7%,  and  aver- 
aging 11.7%.  For  the  Sfax  instances,  the  best  solutions  are  obtained 
with  the  current  fleet  composition,  i.e„  by  using  two  vehicles  of  each 
type.  One  possible  explanation  is  that  compartments  equipped  with 
debit  meters  enable  the  collection  of  small  quantities  of  the  same 
product  and  the  segregation  of  non-mixable  products. 

Finally,  we  compare  our  results  to  the  solution  currently 
applied  by  our  industrial  partner.  Table  4 indicates  that  the 
proposed  method  provides  an  improvement  over  the  current 
solution  designed  manually  by  the  dispatcher  since  it  reduces 
the  overall  costs  for  all  the  data  sets.  Our  solutions  minimize  the 
distance  traveled  with  an  improvement  of  up  to  7%  and  optimize 
the  products  assignment  process  to  avoid  unnecessary  cleaning 
costs.  For  example,  our  solution  for  the  instance  *Kair_27_3_0 
illustrates  a situation  in  which  one  can  save  more  by  performing 
an  extra  tour  and  overcome  the  cost  of  the  cleaning  activity. 
However,  the  same  number  of  vehicles  is  needed  to  cover  all  the 
producers’  locations.  We  have  run  further  tests  with  a hierarchical 
objective  function,  which  first  minimizes  the  required  number  of 
vehicles,  and  then  the  routing  and  cleaning  costs.  These  tests 
reveal  that  the  current  number  of  vehicles  used  by  our  industrial 
partner  is  in  fact  optimal.  Specifically,  comparing  the  structure  of 
the  optimized  solution  with  that  of  the  manual  solution  reveals 
that  the  managers  tend  to  design  tours  starting  by  visiting  the 
producers  which  are  near  to  the  depot,  but  this  intuitive  choice  is 
not  necessarily  optimal.  We  have  also  noticed  that  for  the  periods 
including  few  offers,  the  solutions  found  by  the  managers  were 
optimal  since  there  are  relatively  few  possible  tours.  For  the 
managers,  optimizing  both  the  collection  routes  and  the  fleet 
composition  is  important,  but  difficult  to  achieve  through  manual 
methods.  They  were  asked  to  evaluate  the  solutions  obtained  by 
our  methodology  and  declared  themselves  very  satisfied  with  the 
results. 

In  order  to  test  the  limits  of  the  proposed  algorithm,  we  have 
created  a set  of  15  additional  instances.  These  artificial  instances 
are  generated  similarly  to  real-life  instances  but  on  a larger  scale 
and  include  up  to  19  producers  and  105  requests.  We  have  run  the 
proposed  algorithm  with  all  the  proposed  cuts  (1 )— (13),  (16)— (24) 
over  the  generated  data  sets.  The  computational  experiments 
show  that  instances  including  10  producers  and  60  requests 
cannot  be  solved  to  optimality  within  3 h of  computing  time. 


We  have  tackled  a real-world  and  rich  multi-compartment  f ^ 
vehicle  routing  problem  arising  in  the  olive  oil  collection  industry.  ~ cf 
This  practical  application  concerns  the  pick-up  of  one  or  more  "§Z3 
commodities  from  a set  of  geographically  scattered  producers  in  "§fn 
the  center  of  Tunisia.  We  have  presented  a mathematical  model  cj>m 
including  known  and  new  valid  inequalities,  as  well  as  a branch- 
and-cut  algorithm  for  its  solution.  We  have  shown  how  our  model 
can  be  adapted  to  cover  a variety  of  other  industries  and  assump- 
tions. We  have  tested  the  algorithm  on  several  sets  of  realistic 
instances.  Our  experimental  results  show  the  effectiveness  of  the 
proposed  algorithm.  Extensive  tests  have  enabled  us  to  generate 
solutions  that  can  help  support  decision  making  at  the  tactical 
level,  i.e.,  the  purchase  of  new  vehicles,  and  at  the  operational 
level,  i.e.,  the  redesign  of  vehicle  loading  and  of  collection  routes. 
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